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I.  RESEARCH  OBJECTIVES  AND  APPROACH 


The  basic  objective  of  this  research  program  is  to  obtain  a  clear  phy¬ 
sical  and  theoretical  understanding  of  the  dynamics  of  the  turbulent  boun¬ 
dary  layer  which  will  ultimately  provide  improved  models  for  the  turbulence 
quantities  in  the  time-mean  boundary-! aye r  equations.  The  long  range  goals 
of  the  program  continue  to  be  both  the  improvement  of  the  turbulent  boundary- 
layer  prediction  methods  and  development  of  rational  methods  for  control 
and/or  modification  of  turbulent  boundary  layer  behavior. 

This  project  combines  both  experimental  video  flow  visualization  studies 
and  theoretical  investigations  of  a  series  of  phenomenological  and  theoreti¬ 
cal  models  based  upon  the  three-dimensional  details  of  convected,  coherent 
structural  elements  of  a  turbulent  flow  as  it  interacts  with  a  solid  surface. 
The  experimental  program  is  considering  a  range  of  sub-problems  including 
the  effect  of  surface  modification  on  low  speed  streak  formation  and  drag, 
and  the  effect  of  vortex  loop  interaction  with  a  solid  boundary.  To  augment 
visual  studies,  a  computerized  interface  with  the  video  system  has  been 
developed  which  allows  quantitative  data  to  be  obtained  from  flow  visualiza¬ 
tion  pictures.  The  specific  thrust  of  the  theoretical  studies  has  been 
focussed  on  three  areas:  1)  how  two-  and  three-dimensional  vortex  structures 
interact  with  wall  boundary  layers,  2)  the  development  of  a  new  type  of 
prediction  method  for  two-dimensional  turbulent  boundary-layer  flows,  and 
3)  improvement  in  numerical  techniques  for  solving  parabolic,  boundary- layer 
eq uations . 

II.  STATUS  OF  RESEARCH 


Experimental  Program 

During  the  past  year,  the  Lehigh  experimental  program  on  Coherent  Struc¬ 
ture  of  Turbulent  Boundary  Layers  has  pursued  investigations  in  four  comple¬ 
mentary  areas:  1)  Identification  and  quantification  of  turbulence  structure 
characteristics ,  2)  Control  of  turbulence  structure,  3)  Simulation  of  hypo¬ 
thesized  turbulence  flow  structures,  and  4)  Recreation  of  3-D  motion  in  a 
turbulent  boundary  layer  using  computer  augmented  display  of  video  information. 
These  studies  have  made  heavy  use  of  high-speed  flow  visualization  to  iden¬ 
tify,  characterize,  and  quantify  the  characteristics  of  the  flows  in  question. 
In  addition,  a  substantial  amount  of  parallel  hot-film  anemometry  data  have 
been  taken  for  quantification  of  flow  structure  effects  and  to  assure  com¬ 
parison  with  accepted  turbulence  characteristics. 

Turbulent  Structure  Characteristics 


Our  studies  of  fully  turbulent  flows  have  centered  on  establishing  the 
characteristics  of  the  near-wall  region  of  a  turbulent  boundary  layer  where 
the  most  organized  behavior,  low-speed  streaks,  occurs.  As  pointed  out  last 
year,  the  non-dimensional  spacing  of  the  low-speed  streaks  have  been  shown 
to  be  universally  present  and  to  maintain  a  constant  spanwise  spacing  of 
A+  =  100  over  a  wide  range  of  Reynolds  number  (Figure  1).  In  addition,  we 
have  established  that  the  distribution  of  this  spanwise  spacing  also  appears 
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to  be  universal,  closely  approximating  a  log  normal  distribution  (Figure  2). 
This  latter  piece  of  information  is  of  particular  importance  to  new  predic¬ 
tion  and  turbulence  simulation  efforts  which  require  statistical  flow  struc¬ 
ture  models  as  input. 

An  additional  discovery,  which  has  substantial  implications  with  regard 
to  the  flow  structure  governing  near-wall  behavior,  is  the  observation  of  the 
"persistence"  of  low-speed  streaks,  i.e.  the  tendency  of  streaks  to  perpetuate 
themselves.  This  self  perpetuation,  despite  continual  disruption  by  the 
quasi -periodic  bursting  behavior,  implies  that  the  flow  structure  which  forms 
the  streak  must  also  be  dependent  upon  the  presence  of  the  streaks  for  its 
origin.  Detailed  visual  studies  indicate  that  this  structure  is  a  series  of 
nested  loop-like  vortices  (Figure  3)  generated  as  a  consequence  of  the  "burst¬ 
ing"  of  a  streak.  A  continuing  study  attempting  to  validate  this  hypothesis 
is  described  below  under  simulation  studies. 

Turbulence  Control 

As  a  result  of  our  observations  of  low-speed  streak  "persistence",  we 
have  been  conducting  a  study  to  examine  if  artificial  modifications  of  the 
flow  surface  in  the  straamwise  direction  can  act  as  "sites"  for  streaks  and 
thus  "lock"  them  in  place.  Using  very  fine  fishing  line  as  our  artificial 
streaks,  we  have  established  that  streak  formation  can  be  stabilized  by  the 
fish  line  "sites".  Hot-film  studies  just  completed  indicate  that  the  boundary 
layer  appears  normal  by  all  comparisons,  except  in  the  near-wall  region. 
Spectral  analyses  to  determine  modification  of  bursting  behavior  and  calcula¬ 
tion  of  induced  changes  in  skip  friction  characteristics  are  still  in  progress. 

This  study  of  streamwise  surface  modification  is  very  significant  since 
it  indicates  that  near-wall  turbulent  structure  can  be  controlled  by  passive 
means.  Thus,  we  should  potentially  be  able  to  determine  surface  geometries 
employing  streamwise  modifications  which  can  reduce  surface  drag  by  modifi¬ 
cation  of  the  mechanism  of  momentum  exchange  at  a  surface.  Work  is  continu¬ 
ing  on  this  study,  expanding  our  investigations  to  riblet  surface  geometries 
which  NASA-Langley  has  shown  to  reduce  surface  drag  by  10%. 

Turbulence  Simulation 

The  premise  upon  which  the  Lehigh  turbulence  program  and  a  number  of 
other  programs  across  the  country  are  based  is  that  complex  turbulent  flows 
are  "built"  of  different  "flow  structures."  However,  determining  the  basic 
flow  structures  and  extracting  their  role  in  the  complex  entanglement  of 
fluid  motion  which  is  turbulence  is  an  extremely  difficult  task,  complicated 
by  the  fact  that  we  do  not  understand  the  behavior  of  some  of  the  simpler 
flow  structures  turbulence  has  been  speculated  to  be  comprised  of.  There¬ 
fore,  we  have  undertaken  investigations  of  the  behavior  of  two  simple  flow 
structures  which  appear  to  have  the  potential  for  modeling  some  of  the  key 
characteristics  of  turbulence. 

The  most  mature  of  these  studies  is  of  the  interaction  of  a  ring  vortex 
as  it  impacts  a  solid  surface.  As  pointed  out  last  year,  a  complex  inviscid- 
viscous  interaction  takes  place  which  generates  secondary  and  tertiary  vor¬ 
tices.  The  resulting  group  of  vortices  then  interact  in  a  very  three- 
dimensional,  but  symmetric  fashion  to  rapidly  disperse  the  vorticity  of  the 


initial  vortex.  During  the  past  year,  over  90  different  initial  conditions 
have  been  examined,  indicating  the  complexities  which  can  arise  from  such 
interactions.  Figure  4  shows  a  top  view  of  a  vortex  impact  with  a  surface 
covered  with  dye  which  illustrates  the  extreme  three  dimensionality  of  the 
interaction  process.  Figure  5  shows  one  of  our  dramatic  discoveries,  i.e. 
that  an  impacting  vortex  can  create  a  secondary  vortex  ring  which  "rebounds" 
back  in  opposition  to  the  direction  of  travel  of  the  original  vortex  ring. 

It  is  felt  that  this  effect  may  have  implications,  although  indirectly,  to 
the  bursting  and  ejection  behavior  observed  in  a  turbulent  boundary  layer. 

As  mentioned  above,  our  turbulence  studies  have  suggested  that  nested 
loop-vortices  appear  to  be  generated  by  the  bursting  process  and  to  result  in 
the  formation  of  low-speed  streaks.  To  examine  this  hypothesis,  we  have 
developed  a  technique  for  generating  very  consistent  loop  vortices  similar 
in  appearance  to  the  schematic  shown  in  Figure  3.  It  was  found  that  under 
the  proper  conditions,  loop  vortices  will  be  shed  in  a  very  periodic  fashion 
from  a  hemisphere  placed  on  a  surface  beneath  a  laminar  boundary  layer.  Very 
extensive  visualization  and  hot-film  probe  studies  have  been  done  of  these 
loops,  and  the  results  have  proven  extremely  encouraging.  Thus  far  it  has 
been  determined  that  many  of  the  visual  patterns  observed  in  the  near-wall 
region  of  a  turbulent  boundary  layer  can  be  observed  in  the  wake  of  the 
hemisphere. 

This  study  of  loop  vortices  is  continuing  in  order  to  better  establish 
the  parallel  between  these  systematically  generated  structures  in  a  laminar 
boundary  layer  and  the  less  well  defined  structures  arising  from  the  near- 
wall  of  a  fully  turbulent  boundary  layer. 

Recreation  of  3-0  Motion 


In  our  previous  work,  the  dual -camera  video  system  was  used  in  conjunc¬ 
tion  with  a  fiber-optic  lens  •  obtain  combined  top-end  views  of  single  real¬ 
ization  hydrogen  bubble-lines  in  the  near-wall  region  of  a  turbulent  boundary 
layer.  From  these  combined  views,  we  were  able  to  establish  a  number  of 
characteristics  regarding  behavior  in  the  near-wall  region,  particularly  the 
characteristics  and  potential  formation  mechanism  for  low-speed  streaks. 
However,  the  presentation  of  the  two-view,  time-sequence  behavior  was  awkward 
(requiring  numerous  pictures).  Additionally,  the  two-views  could  not  be 
taken  in  true  orthogonal  perspective,  which  provoked  some  confusion  for  the 
unfamiliar  observer. 

To  rectify  these  shortcomings,  and  to  provide  the  potential  for  reduc¬ 
tion  of  these  two-view  sequences  to  quantitative  velocity  field  information, 
a  system  for  digitizing  the  bubble-line  sequences  into  a  computer-aided 
display  system  has  been  developed.  Using  this  system,  the  Lagrangian  motion 
of  a  single  */Ubble-line  in  time  and  space  can  be  displayed  in  proper  ortho¬ 
graphic  perspective  (computer  corrected  for  viewing  orientation)  as  shown  in 
Figure  6.  Since  the  digitized  bubble-line  information  is  stored  as  a  three- 
dimensional  matrix,  any  view  of  the  recreated  bubble-line  motion  and  deforma¬ 
tion  can  also  be  displayed,  as  shown  by  the  oblique  view  in  Figure  6.  In  addi¬ 
tion,  by  computer  manipulation,  orthogonal  splines  can  be  matched  to  the 
original  bubble-lines  to  create  a  grid-like  "surface"  of  the  bubble-line 
motion.  Such  a  surface,  created  from  the  bubble-line  data  of  Figure  6,  is 
shown  in  Figure  7. 
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Presently  this  computer  display  process  is  being  used  to  examine  a 
number  of  near-wall  sequences.  In  addition,  modifications  to  both  the 
display  program  and  to  the  visualization  system  are  being  done  which  will 
allow  velocity  field  information  to  be  derived  from  the  same  input  as  was 
used  to  create  Figures  6  and  7. 


ANALYTICAL  PROGRAM 


The  objective  of  this  report  is  to  review  the  progress  made  over  the 
past  year  as  well  as  to  delineate  the  current  direction  of  the  program 
and  the  expected  results  that  will  be  obtained  by  the  end  of  the  current 
contract  year  (April,  1983).  The  analytical  program  is  directed  in  three 
separate  but  related  phases.  The  first  of  these  areas  is  development  of 
improved  numerical  methods  for  parabolic  partial  differential  equations. 
This  area  is  an  outgrowth  of  previous  work  in  the  development  of  turbulent 
boundary  layer  prediction  methods  where  the  need  to  improve  the  efficiency 
of  the  calculation  procedures  became  apparent.  The  overall  goal  of  this 
portion  of  the  program  is  to  develop  a  numerical  procedure  for  nonlinear 
pa»"»holic  equations  which  is  both  fourth  accurate  and  allows  an  arbitrary 
non-uniform  mesh  spacing.  This  latter  requirement  is  particularly  impor¬ 
tant  in  the  calculation  of  turbulent  boundary  layer  flows  where  packing 
of  the  mesh  near  the  wall  is  necessary  to  achieve  good  accuracy.  A  fourth 
order  method  is  important  to  substantially  reduce  computer  storage  require¬ 
ments  as  well  as  computation  times.  At  this  stage,  it  appears  the  goal 
will  be  realized  by  April,  1983. 

During  the  past  year,  Mr.  W.-C.  Lee  completed  work  on  his  Master's 
degree  program  and  his  thesis  work  concerned  the  development  of  a  new 
second  order  accurate  parabolic  equation  method.  The  thesis  is  currently 
in  preparation  as  an  AFOSR  technical  report.  Some  typical  results  are 
illustrated  in  Figure  8  where  the  root-mean-square  error  for  two  new 
calculation  methods  is  compared  with  that  associated  with  two  existing 
second  order  methods  (the  Keller  Box  method  and  the  Crank-Nicolson  method). 
The  example  problem  for  which  the  error  is  illustrated  in  Figure  1  is 
the  Howarth  laminar  boundary  layer  problem*,  this  flow  is  a  boundary  layer 
developing  in  an  adverse  pressure  gradient,  similar  to  a  diffuser- type 
flow.  The  boundary  layer  originates  at  £=0  and  a  separation  point  is  pre¬ 
dicted  at  c  =  0.90;  the  RMS  error  plotted  in  Figure  8  is  the  average  error 
incurred  by  each  method  at  each  5  station  in  a  boundary  layer  integration 
initiated  at  s  =  0.  It  may  be  observed  that  the  new  methods  offer  a  good 
increase  in  accuracy.  The  study  of  Mr.  Lee  is  also  important  because  it 
illustrated  that  spatial  differencing  techniques  developed  for  ordinary 
differential  equations  may  be  carried  over  with  some  additional  ingenuity 
to  parabolic  partial  differential  equations.  During  the  past  year  a 
fourth  order  scheme  for  both  linear  and  nonlinear  equations  has  been 
developed;  this  method  is  restricted  to  uniform  mesh  spacings.  The  results 
of  this  study  are  part  of  the  Ph.D.  thesis  work  of  Mr.  E.A.  Bogucz  and 
have  been  submitted  for  publication.  The  results  of  this  study  are  very 
encouraging  and  illustrate  how  the  number  of  mesh  points  may  be  reduced  by 
two  orders  of  magnitude.  At  present,  there  are  three  remaining  tasks  to 
complete  this  portion  of  the  study;  these  are  development  of: 

(a)  a  fourth  order  method  for  parabolic  partial  differential 
equations  using  a  uniform  spatial  mesh; 

(b)  a  fourth  order  method  for  ordinary  differential  equations 
with  an  arbitrary  non-uniform  mesh; 

(c)  a  fourth  order  method  for  partial  differential  equations 
with  an  arbitrary  non-uniform  mesh. 
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Currently  the  feasibility  of  tasks  (a)  and  (b)  has  been  demonstrated  and 
the  methods  under  development  are  being  tested.  It  is  evident  that  the 
methods  being  produced  here  will  have  wide  applicability  in  other  areas 
of  engineering  and  the  physical  sciences. 

The  second  area  of  present  effort  is  concerned  with  the  development 
of  improved  turbulence  models.  In  previous  years  of  this  contract,  a 
model  for  the  inner  layer  of  the  turbulent  boundary  layer  was  developed 
which  was  based  on  the  observed  coherent  nature  of  the  wall  layer  flow. 
Current  interest  surrounds  development  of  a  model  for  the  outer  layer  and 
in  this  portion  of  the  contract  work,  consideration  is  given  to  obtaining 
a  simple  eddy  viscosity  model.  During  the  past  year,  Mr.  L.J.  Yuhas 
completed  a  Master's  thesis  (which  has  been  submitted  as  an  AFOSR  technical 
report);  in  this  study,  a  procedure  is  described  wherein  the  basic  para¬ 
meters  in  the  eddy  viscosity  model  may  be  correlated  to  account  for  various 
physical  effects;  the  correlations  are  obtained  by  comparing  a  turbulent 
velocity  profile  directly  with  measured  experimental  profile  data.  An 
optimization  computer  code  has  been  developed  to  carry  this  procedure  out 
and  the  method  has  been  applied  to  flows  with  pressure  gradient  and  flows 
with  mainstream  turbulence.  As  an  illustration  of  how  well  the  velocity 
profiles  match,  some  data  (taken  under  an  AFOSR  program  at  United  Techno¬ 
logies  Research  Center,  East  Hartford,  Conn.)  for  mainstream  turbulence 
levels  ranging  from  3.5%  to  6.5%  is  illustrated  in  Figure  9  for  a  zero 
pressure  gradient  flow.  Such  mainstream  turbulence  levels  are  typical  of 
boundary  layer  flows  in  the  gas  turbine  environment.  It  may  be  observed 
in  Figure  9  that  the  profiles  calculated  in  the  optimization  procedure 
represent  the  measured  data  very  well.  From  these  optimizations  it  proved 
possible  to  correlate  a  single  parameter  K  in  the  simple  eddy  viscosity 
formula.  The  quadratic  correlation  is  illustrated  in  Figure  TO  and  may 
now  be  used  in  the  eddy  viscosity  formula  in  a  prediction  method. 

A  second  approach  to  the  problem  of  model  development  of  the  outer 
layer  of  the  turbulent  boundary  layer  constitutes  the  third  and  a  major 
phase  of  the  research.  In  this  portion  of  the  contract,  physical  mechanisms 
and  the  underlying  causes  of  the  observed  structure  and  bursting  event  in 
turbulent  boundary  layers  are  under  investigation.  During  the  past  year, 
a  study  has  been  completed  on  the  viscous  effects  that  are  due  to  be 
expected  due  to  a  pair  of  counter-rotating  vortices  above  a  plane  wall. 

Such  vortices  are  observed  in  the  wall  region  of  a  turbulent  boundary  layer 
and  are  believed  by  some  authors  to  play  an  important  role  in  the  dynamics 
of  the  production  of  turbulence.  Moreover  the  longitudinal-1-  vortices  in  a 
turbulent  boundary  layer  are  observed  to  persist  for  long  periods  of  time 
and  apparently  do  not  migrate  outwards  as  a  purely  inviscid  theory  would 
suggest  they  should.  Consequently  it  was  decided  to  undertake  a  study  to 
elucidate  the  viscous  effects  induced  near  a  wall  by  a  pair  of  counter¬ 
rotating  vortices.  The  study  that  was  completed  this  year  was  for  a  pair 
of  two-dimensional  vortices  which  initially  are  at  the  same  height  above 
a  plane  wall;  inviscid  theory  predicts  that  both  vortices  will  move  on  the 


+  In  fact,  the  observed  vortices  cannot  be  entirely  longitudinal  and  must 
be  portions  of  an  elongated  three  dimensional  loop. 
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hyperbolic  illustrated  in  Figure  11.  The  vortices  will  either  move  away 
from  the  wall  or  toward  it.  depending  on  the  sense  of  rotation  of  the  vor¬ 
tices.  The  flow  is  symmetric  about  a  line  bisecting  the  vortex  paths.  In 
all  cases  considered,  a  boundary  layer  separation  occurs  leading  to  the 
creation  of  a  counter-rotating  pair  of  secondary  vortices.  In  the  case 
of  an  upward  moving  pair  of  vortices,  separation  occurs  on  the  wall  inboard 
of  the  rising  vortices.  A  typical  case  is  illustrated  in  Figure  12  where 
a  developing  separation  occurs  on  the  symmetry  plane;  note  that  the  illus¬ 
trations  in  Figure  12  are  for  one  half  of  a  symmetric  boundary  layer  flow 
and  that  there  is  a  mirror-image  boundary  layer  flow  for  x<0.  If  the  vor¬ 
tices  are  started  further  apart,  a  separation  occurs  which  again  is  inboard 
of  the  rising  vortices  but  which  is  not  connected  with  the  mirror  image 
on  the  symmetry  plane.  Such  a  case  is  illustrated  in  Figure  13.  Note 
that  in  Figures  5  and  6  the  bold  arrows  at  the  top  of  the  figures  denote 
the  current  streanwise  location  of  the  parent  vortex  in  the  inviscid  flow; 
the  smaller  arrow  indicates  the  initial  starting  location  of  the  vortex. 

In  the  case  of  the  upward  moving  pair,  the  eruption  of  the  secondary  vor¬ 
tices  from  the  boundary  layer  will  arrest  the  upward  and  inward  movement 
of  the  parent  vortices.  For  a  downward  moving  pair,  it  was  determined  that 
the  creation  of  secondary  vortices  in  the  boundary  layer  also  occurs  but 
now  outboard  of  the  parent  vortices.  Consequently  in  this  case  as  well, 
the  eruption  of  the  secondary  vortices  will  act  to  arrest  the  motion  of 
the  parents.  This  study  is  part  of  the  Ph.D.  program  of  Hr.  S.  Ersoy  and 
a  paper  describing  the  complete  results  is  currently  in  preparation. 

At  present,  the  vortex  interaction  part  of  the  program  is  directed 
towards  carrying  out  the  following  major  tasks: 

(a)  the  calculation  of  the  boundary  layer  induced  by  a  pair  of 
vortices  moving  toward  a  wall  at  an  angle 

(b)  the  calculation  of  the  inviscid  flow  due  to  a  three  dimensional 
loop  vortex  and  the  boundary  layer  flow  induced  by  the  motion 
of  the  loop. 

At  this  time,  work  is  well  along  the  way  for  task  (a).  In  Figure  14,  a 
typical  trajectory  for  a  vortex  pair  approaching  a  wall  obliquely  is  given. 
The  boundary  layer  flow  induced  by  an  asymmetric  pair  of  vortex  distur¬ 
bances  is  much  more  complex  than  the  previous  problems  considered  in  this 
contract.  For  this  reason,  it  has  required  a  considerable  effort  to 
deduce  the  appropriate  analytic  transformations  and  numerical  methods 
required  to  handle  the  problem.  Some  initial  results  have  been  produced 
and  in  Figure  15  the  boundary  layer  flow  due  to  the  pair  in  Figure  14  is 
plotted  at  a  certain  stage  in  the  boundary  layer  development.  Note  that 
a  secondary  eddy  has  been  created  in  the  region  near  x=-3  by  the  lower  of 
the  parent  vortices  (in  Figure  14).  There  is  a  tendancy  for  a  weaker  secon¬ 
dary  separation  near  x=2.5  but  at  the  stage  of  development  in  Figure  15, 
the  second  secondary  separation  has  not  yet  occurred.  Over  the  next  year 
a  number  of  such  cases  will  be  considered. 

Work  has  started  on  task  (b)  and  at  this  stage  is  focused  primarily  on 
the  calculation  of  the  inviscid  three  dimensional  loop  motion.  The  boundary 
layer  flow  in  this  case  is  somewhat  more  complex  than  case  (a);  however, 
many  of  the  techniques  of  solution  developed  for  task  (a)  can  and  will  be 
adapted  to  the  vortex  loop  problem  in  the  next  contract  year. 
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Figure  2.  Probability  density  histograms  of  spanwise  streak  spacing  at  y  =5. 

a)  Re  a  1490,  n  =  437;  b)  Re„  8  5830,  n  8  411.  - — ,  lognormal  prob¬ 
ability  density  distribution  for  corresponding  A  and  . 


Figure  5  Secondary  vortex  ejection  (side  view,  dye  placed 
on  surface).  Vortex  ID,  #52.265,  Reo=3000.  Sequence  continued 
on  next  page. 


bubble-line  motion. 


Crank-Nicolson  Method 


Figure  8.  Comparison  of  RMS  error  for  the  Howarth  flow  problem 


Figure  9.  Velocity  profile  comparisons  for  two  parameter 
optimization  on  S  and  K  {with  <=0.447B9)  for 
mainstream  turbulence  flaw. 


K  parameter  correlation  obtained  from  one  parameter  fits  of  the  composite 
similarity  profile  to  constant  pressure  test  data  with  mainstream  turbulence 


moving  vortex  pair  started  at 
Xo=Vo=*/2  at  t=0.65 . 


